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Abstract 

In this paper, we extend Walsh's stochastic integral with respect to 
a Gaussian noise, white in time and with some homogeneous spatial cor- 
relation, in order to be able to integrate some random measure-valued 
processes. This extension turns out to be equivalent to Dalang's one. 
Then we study existence and regularity of the density of the probability 
law for the real- valued mild solution to a general second order stochastic 
partial differential equation driven by such a noise. For this, we apply the 
techniques of the Malliavin calculus. Our results apply to the case of the 
stochastic heat equation in any space dimension and the stochastic wave 
equation in space dimension d = 1,2,3. Moreover, for these particular 
examples, known results in the literature have been improved. 
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1 Introduction 



The purpose of this paper is to study the probabiUty law of the real-valued 
solution of the following general class of stochastic partial differential equations: 

Lu{t, x) = a{u{t, x))W{t, x) + b{u{t, x)), (1.1) 

t > 0, a; S M'^, where L denotes a second order differential operator, and we 
impose the initial conditions 

u{Q,x)^^{Q,x)^Q. (1.2) 

The coefficients a and h are some real- valued functions and W{t^ x) is the formal 
notation for some Gaussian random perturbation defined on some probability 
space. We will assume that it is white in time and with a homogeneous spatial 
correlation given by a function /, and we denote by (^t)t the filtration gener- 
ated by W (see Section [2] for the precise definition of this noise). 

By definition, the solution to Equation ()l.ip is an adapted stochastic process 

{u{t,x), {t,x) e [0,r] X R''} satisfying 

u{t,x)^ / ^{t ~ s,x - y)a{u{s,y))W{ds,dy) 
Jo Js.'' 

+ 11 b{u{t~s,x~y))T{s,dy)ds, (1.3) 

where F denotes the fundamental solution associated to Lu ~ 0. As it is going 
to be clarified later on, for all t e [0, T], we suppose that r{t) is a non-negative 
measure on (see Hypothesis A). The stochastic integral appearing in formula 
(|1.3p requires some care because the integrand is a measure. For integrands that 
are real- valued functions, that is {t, x) i-^ F(i, x) e R, the stochastic integral was 
defined by Walsh in [16]. Then, in order to deal with SPDEs whose associated 
fundamental solution is a generalised function, Dalang ^ extended Walsh's 
stochastic integral and covered, for instance, the case of the wave equation in 
dimension three. 

However, in the first part of this paper we will give, in a general setting, a 
definition of the stochastic integral with respect to W using the techniques of the 
stochastic integration with respect to a cylindrical Brownian motion (see, for in- 
stance, [5]). As we will see, this integral will turn out to be equivalent to Dalang's 
extension ([3j) when the integrand is of the form G := T{t — •, a; — *)Z{-, *), for 
certain stochastic processes Z. More precisely, we will show that random ele- 
ments of this latter form may be integrated with respect to W using a localising 
procedure: first we will assume that Z has bounded trajectories and then we 
will identify G as the weak limit of some sequence (FZjv)Ar, where any have 
bounded paths, almost surely (see Lemma [3?2l and Proposition 13. Sp . 
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We should mention at this point that solutions to stochastic partial differ- 
ential equations of the form (|1.1[) have been largely studied during the last two 
decades. For instance, for the case of the stochastic heat and wave equations, 



The second part of the paper is devoted to study the probability law of the 
random variable u{t,x), for any fixed {t,x) G (0,r] x R'^. This will be done 
using the techniques of the so-called Malliavin calculus. The aim is two-fold: 

First, we prove that u{t, x) has an absolutely continuous law with respect to 
Lebesgue measure on R, provided that, among other assumptions, the differen- 
tial operator L and the spatial correlation / are related as follows: 



where /z is a non-negative tempered measure such that its Fourier transform is 
/ (see Theorem 15.21 for the precise statement). This result provides a generali- 
sation of Theorem 3 in [13] , where the authors deal with the three-dimensional 
stochastic wave equation and a slightly stronger condition than (|1.4p is assumed. 
In order to prove that the law ofu{t, x) has a density, we apply Bouleau-Hirsch's 
criterion. Indeed, to prove the Malliavin regularity of the solution, we take ad- 
vantage of the results in [13j and we fully identify the initial condition of the 
stochastic equation satisfied by the Malliavin derivative of u{t,x). This is an 
important point in order to show that the Malliavin matrix is invertible almost 
surely. Eventually, we point out that (II. 4|) is also a sufficient condition to have 
the stochastic integral in (|1.3|1 well defined. 

Secondly, we prove that the law of u{t,x) has an infinitely differentiable 
density with respect to Lebesgue measure on R (see Theorem I6.2|) . To obtain 
this result, we show that u{t, x) is infinitely differentiable in the Malliavin sense 
and that the inverse of the Malliavin matrix has moments of all order. For the 
latter to be achieved, we need to impose a lower bound of the integral in (|1.4p 
in terms of a certain power of T (see (|6.29p ). We should mention that Theorem 
16.21 in Section [6] provides an improvement of Theorem 3 in [14] for the case of 
the three-dimensional stochastic wave equation, since in this latter reference the 
integrability conditions concerning the Fourier transform of F are much more 
involved (see also !T5|)). Moreover, it is worth mentioning that Theorem 16.21 
also generalises known results on existence and smoothness of the density for 
the case of the stochastic heat and wave equation with dimensions d > 1 and 
d= 1,2, respectively (see [HHHTHH]). 

The paper is organised as follows. In the next Section |21 we present some 
preliminaries concerning the random perturbation in Equation (jl.ip as well as 
the main hypothesis on the fundamental solution F and the space correlation 
/. We extend Walsh's stochastic integral in order to cover the case of measure- 
valued integrands in Section |3l we also caracterise some measure- valued random 



we refer to [M HH li 13 HIl [IH] ■ 




(1.4) 
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elements that can be integrated with respect to W and we sketch the construc- 
tion of the integral in a Hilbert-valued setting. In Section |4l a theorem on 
existence and uniqueness of solution for Equation (jl.ip is stated; we also deal 
with some particular examples of differential operators L, namely the heat and 
wave equations. Sections [5] and [S] are devoted, respectively, to the existence 
and regularity of the density for the probability law of the solution to p.ip . At 
the very beginning of Section [5l we introduce the main tools of the Malliavin 
calculus needed along the paper (we refer to [5] for a complete account on the 
topic). 

Throughout the paper we use the notation C for any positive real constant, 
independently of its value. 



2 Preliminaries 

Recall that we are interested in the following general class of stochastic partial 
differential equations: 

Lu{t, x) = a{u{t, x))W{t, x) + h{u{t, x)), 

with i > 0, x e W^, L denotes a second order differential operator, and we 
consider vanishing initial conditions (|1.2p . The Gaussian random perturbation 
is described as follows: is a zero mean Gaussian family of random variables 
{W{ip), (fi e C,f=(R'*+i)}, defined in a complete probability space (O, J^, P), with 
covariance 

E{W{^)W{ij))^ I I I ^(t,x)f{x-y)^{t,y)dxdydt, (2.5) 
Jo Jr'^ Jr'' 

where / is a non-negative continuous function of M.'^ \ {0} such that it is the 
Fourier transform of a non- negative definite tempered measure /i on M'*. That 
is, 



fix) = / exp(-27ri x ■ Ofi{d^) 
and there is an integer m > 1 such that 

(1 + ieP)~™A*(dO < 



Then, the covariance (12. 5p can also be written, using Fourier transform, as 



E{W{ip)W{4>)) - / / Tipit)iOTmiOKdOdt. 
Jo Jr'' 

The main assumption on the differential operator L may be summarised as fol- 
lows: 
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Hypothesis A. The fundamental solution to Lu = 0, denoted by F, is a non- 
negative measure of the form T{t, dx)dt such that for all T > 

sup T{t,W^) <Ct<oo (2.6) 

0<t<T 

and ^ 

/ / \Trm)fKdOdt<oo. (2.7) 

Jo JR<* 

The completion of the Schwartz space iS(M'') of rapidly decreasing C°° func- 
tions, endowed with the inner product 

{cp,i;)H= f f ip{x)f{x-y)ip{y)dxdy= f Tip{OWiO^^{dO, 

<p,ip €z S(R'^), is denoted by H. Notice that Ti may contain distributions. Set 
HT = i'([0,T];W). 

3 Stochastic integrals 

In this section Walsh's stochastic integral with respect to martingale measures 
will be extended to more general integrands, namely the class of square inte- 
grable 7i-valued predictable processes. The extension will be performed in the 
infinite dimensional setting described by Da Prato and Zabczyk in [5] . Then, we 
will give non-trivial examples of integrands, which will be some measure- valued 
random elements. We will briefly recall the extension of the stochastic integral 
in a Hilbert-valued setting. This will be needed to give a rigorous meaning to 
the stochastic evolution equations satisfied by the Malliavin derivatives of the 
solution of (II. 3p . 

Extension of Walsh's stochastic integral 

Fix a time interval [0, T]. The Gaussian family {W{ip), ip £ C^([0, T] x R"*)} 
can be extended to the completion Ht = ^^([0, T];H) of the space C^([0, T] x 
R"*) under the scalar product 

Jo Jw 

We will also denote by W{g) the Gaussian random variable associated with an 
element g € L^{[0,T];n). 

Set Wtih) = Wil[o,t]h) for any i > and h e H. Then, {Wt,t e [0,T]} is 
a cylindrical Wiener process in the Hilbert space H. That is, for any h ^ H, 
{Wt{h),t G [0, r]} is a Brownian motion with variance and 

E{Wt{h)W,{g) = isAt) {h,g)^. 
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Let Tt be the cr-field generated by the random variables {Wsiji),h G H,0 < 
s < t} and the P-nuU sets. We define the predictable cr-field as the cr-ficld in 
n X [0, T] generated by the sets {{s,t] x A,0 < s < t < T, A e T^}. 

Then (see, for instance, [S]), we can define the stochastic integral of H- 
valued square integrable predictable processes. For any predictable process g G 
L^{n X [0,T]]H) we denote its integral with respect to the cylindrical Wiener 
process W by 

/ [ gdW = g-W, (3.8) 

Jo JR-i 

and we have the isometry property 




Remark 3.1 Under the standing assumptions, using an approximation proce- 
dure by means of test functions, one proves that the space Ti contains the indi- 
cator functions of bounded Borel sets (for details see ^ or p. 13). Then, 
Mt{A) := W^(1[Q defines a martingale measure associated to the noise W in 
the sense of Walsh (see fTBjj and fS^) and the stochastic integral iS. 8\) coincides 
with the integral defined in the work of Dalang f^. 

Example of integrands 

We aim now to provide useful examples of random distributions which belong 
to the space L^(ilx [0, T];Ti.). Before stating the result, we consider the following 
lemma: 

Lemma 3.2 Assume that T satisfies Hypothesis A. Let g be a bounded Borel 
function on [0,T] x M.'^. Then gT 6 Ht, and 

\\9nn^<\\9\\l r f immoi' t^idodt. 

Jo Jr'' 

Proof. We can decompose g into the difference g+ — g~ of two nonnegative 
bounded Borel functions. Thus, without any loss of generality we can assume 
that g is nonnegative. Moreover, we observe that gT also satisfies conditions 
(|2.6p and (|2.7p . Indeed, to prove the latter condition, we consider an approx- 
imation of the identity {ipn)n defined as follows: let ip S C^{M.'^) such that 
V' > 0, the support of i/' is contained in the unit ball of K'' and ip{x)dx = 1. 
Set r„(t) = r(i) and J„(t) * {gT{t)). Then, for all t e [0,T], r„(t) 

and Jnit) belong to S{W^) C U, and |:^r„(t)| < \TV{t)\. Besides, since F is 
non-negative, we have that Jn{t){0 < llff ||oorn(^)(0> fo^' ^-i^y C G Thus, by 
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Fatou's lemma 



<liminf/ / / Jn{t,x)f{x - y)Jnit,y)dxdydt 

< llffll^liminf / / |^r„(i)(0|V(de)rfi 

< Ml [ [ \mt)iofKdOdt<^. 

Jo Jr-^ 



The fact that gT satisfies conditions (|2.6p and (|2.7p let us reduce the proof to 
the case where g = 1- 

We consider the regularisation (r„)„ of F defined above. Condition ()2.7p 
implies that F„(t) belongs to Tix and it has a uniformly bounded norm, so it 
converges weakly to some element h G Ht- We claim that h = T, and this is a 
consequence of the fact that, owing to the definition of F„(t), for any ip e 5(]R''), 
and for any < s < t < T we have 



{h,ip)^dr ^ lim / {Tn(r),ip)^dr = I {T{r),ip)^dr. 

J s J s 

More precisely, it holds that 

{rn{r),(p)^dr ^ J r{r,dz)^J iljn{x - z)F{x)dx^ dr, 

where F{x) :— J^^a f{x — y)ip{y)dy. Observe that the hypothesis on / and (p 
imply that F is continuous and lim^j^oo F{x) = 0. Hence, it turns out that we 
can apply the Bounded Convergence Theorem, so that we end up with 

lim / (F„(r),v?)„dr = / / / T{r,dz)f{z-y)(p{y)dydr, 

and this let us identify h with F in TCt- □ 

This lemma allows us to prove the following result, which give examples of 
random distributions that can be integrated with respect to W. 

Proposition 3.3 Assume thatT satisfies Hypothesis A. Let Z — {Z{t,x), {t,x) £ 
[0, T] X R'^} be a predictable process such that 

Cz sup E{\Z{t, x)\^) <oo. 
{t,x)e[a,T] xW 

Then, the random element G — G(t,dx) = Z{t, x)T{t, dx) is a predictable pro- 
cess in the space L^{fl x [0,T];7i). 
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Proof. For any > 1 define 

ZN{t,x) = Z{t,x)li\z(t,x)\<N}- 

Clearly, Z^r is a predictable process with bounded trajectories. Thus, by the 
previous lemma, GN{t, x) := ZN{t, x)T{t, dx) is a predictable process in L'^{n x 
[Q,T]]'H). Let {JN.n{t))n C be the regularisation of G'Ar(i) by means 

of an approximation of the identity {ipn)n, as it has been defined in the proof 
of Lemma [3.21 This let us prove that the norm of Gat in L^(ri x [0,T];H) is 
uniformly bounded because 

<liminf£' / / / JN,nit,x)fix - y)JN.n{t,y)dxdydt] 
< Czliminfi; ( / / / Tn{t,x)f{x-y)Tn{t,y)dxdydt\ 



<Cz [ I \TT{tm? tJ^id^dt < oo. 

Jo JR<* 



Recall that r„(i) = * r(i). Therefore, Gat converges weakly to a predictable 
process G in L'^{n x [0, T];H). We claim that G = G. In fact, for any ip e S{R'^), 
for any < s < t < T and for any B g JF^, we can argue similarly as in the very 
last part of the proof of Lemma 13.21 and obtain 



E(lsf{Gir),^)^dr) 



liia^E yis J {GN{r),(p)-^dr ^ 
lim E[1b / / T{r,dx)Z{r, 

x)lnz(rt x)\<N} f {x — y)f{y)dxdydr J 
E[1b / / r{r,dx)Z{r,x)f{x - y)if{y)dxdydr 



so we can identify G with G. □ 

Remark 3.4 As a consequence of Proposition \3.3\ we are able to define the 
stochastic integral of G = ZT with respect to W: 

G-W = r I G{s,y)W{ds,dy)^ C ( V{s,y)Z{s,y)W{ds,dy). 



JVL-i 
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Moreover, using the same ideas as in \18f (see also Theorem 1.2.5), one 

can obtain hounds for the U\Q) — norm of G ■ W . More precisely, suppose that 



for some p > 2. Then 



E{\G-W\P) = E 



sup E{\Z{t, x)\P) < oo, 

(t,a:)e[0,T]xK<' 



G{s,y)W{ds,dy) 



wher 



<C,[utY^-' r (snvE{\Z{s,x)\P)] j \TT{sm\' ^,{di)ds, (3.9) 

\rr{tm?ix{d£,)dt. 



T 

Jm.'^ 



Remark 3.5 Since the noise's correlation is of the form {x,y) i— > f{x — y), it 
is natural to be interested in spatially homogeneous situations. Indeed, suppose 
that we add the following hypothesis on the process Z: for all s € [0,T] and 
x,y G M.'^ we have 

E[Z{s, x)Z{s, y)) ^ E{Z{s, 0)Z(,s, y - x)). 

Then, owing to f^, p. 10, we may construct a non-negative tempered measure 
/if on such that 

l|G|U.(0;«.) = E{\G ■ W\^) ^ f I \TT{sm?l^'{dOds. 

Jo Jw^ 



As we will see in the next section, the main examples of deterministic mea- 
sures r will correspond to fundamental solutions associated to second order 
differential operators. First, we sketch the construction of the stochastic inte- 
gral in a Hilbert- valued setting, that is when the process Z takes values in some 
Hilbert space, usually different from H. 

Hilbert-valued stochastic integrals 

Let be a separable real Hilbert space with inner-product and norm denoted 
by {■,-)a and || • |U, respectively. Let K = {K{t,x),{t,x) G [0, T] x R''} be an 
^—valued predictable process satisfying the following condition: 

sup E{\\K{t,x)\f_^) <oo. (3.10) 

{t,x)<£[0,T]xR'' 

Our purpose is to define the stochastic integral of elements of the form TK = 
r(i, dx)K{t, x) eL^iflx [0, T]-n(g) A). 
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Let (ej , j > 0) be a complete orthonormal system of A. Set {t, x) = 
{K{t,x),ej)A, (t,x) e [0,T] X W^. According to Proposition [3^ for any j > 
the element = G^{t,x) = T{t,dx)Ki {t,x) belongs to L'^{n x [0,T];W) and, 
therefore, we may integrate it with respect to the noise W: 

G'-W^ [ ( T{s,y)K^{s,y)W{ds,dy). 

Jo JR'' 

We define, for G = TK, 

G W ■.= Y,G^ -V^- 

Owing to (|3.10p and Proposition 13.31 it can be proved that the above series is 
convergent and therefore G ■ W defines an element of L'^{n;A) (see also [13], 
Remark 1). Moreover, using the same arguments as for the proof of p.9p . we 
have the following bound for the moments of G ■ in ^: 

E{\\G-W\f^) 

<C^{vtY--' I sup E{\\K{s,x)\f_^) [ \Tr{sm\^fi{d^)ds, (3.11) 

Jo x£R'i JR'^ 

for all p > 2. 

4 Existence and uniqueness of solutions 

Recall that a solution to Equation (jl.ip is a real-valued adapted stochastic 
process u = {u{t,x), {t,x) G [0,T] x R"*} satisfying 

u{t,x)= / r{t - s,x - y)a{u{s,y))W{ds,dy) 
Jo JR'' 

+ [ [ b{u{t - s,x - y))T{s,dy)ds. (4.12) 

Jo JR<* 

We assume that Z = Z{s,y) = a{u{s,y)) satisfy the hypothesis of Proposition 
13.31 so that the stochastic integral on the right hand-side is well-defined. 

We suppose that cr and b are real- valued Lipschitz functions. Under these 
conditions, we may state an existence and uniqueness of solution's theorem: 

Theorem 4.1 Suppose that the fundamental solution F of Lu = satisfies 
Hypothesis A. Then, Equation {^.12^ has a unique solution {u{t,x)^{t,x) G 
[0, r] X W^} which is continuous inL? and satisfies 

sup E{\u{t,x)\P) < oo, 

(t,a:)G[0,T]xR'' 

for all T > and p > 1. 
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For the proof we refer to [3], Theorem 13, where the Walsh-Dalang equivalent 
setting is used. 

Let us now enumerate some examples of differential operators whose associ- 
ated fundamental solution fulfills the hypothesis of Theorem 14.11 

Example 4.2 The wave equation. Let Td be the fundamental solution of the 
wave equation in R"* , that is, Td is the solution of 

ATd = 0, 



with vanishing initial conditions. It is known that for d — 1,2,3, Td is given, 
respectively, by 



'1/2 



1 

ri(^) = 2'^i\x\<t}, 
r,it) = cit^-\xr) 
m) = ^a., 

where at denotes the .surface measure on the three-dimensional sphere of radius 
t. In particular, for each t, Td{t) has compact support. It is important to remark 
that only in these cases Td defines a non-negative measure. Furthermore, for 
all dimensions d>l, we have a unified expression for the Fourier transform of 
Td{t): 

sin(2^t|e|) 
^^'^"^^^^ ^ 2vr|e| ■ 

Elementary estimates show that there are positive constants ci and C2 depending 
on T such that 

c, r-sin^^^^ C2 



i + iep-Jo 47r2|cp -1 + iei 

Therefore, Td satisfies condition \2. 7| ) if and only if 

fiidO 



Example 4.3 The heat equation. LetT be the fundamental solution of the heat 
equation in and with vanishing initial conditions, that is 

dt 2 

Then, T is given by the Gaussian density: 

T{t,x) = {2nt)-''/'expU-^ 
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and 



Because 



we conclude that condition ^2.7^ holds if and only if l{4.13\ l holds. 



Let us express condition (I4.13P in terms of the covariance function /. Indeed, 
as it is pointed out in [3J, condition (I4.13P is always true when d — I; for d = 2, 
(|4T3ll holds if and only if 



5 Existence of density 

In this section we aim to prove that the solution to Equation (|1.3p . at any point 
{t, x) £ (0, T] X R'^, is a random variable whose law admits a density with respect 
to Lebesgue measure on R. For this, we will make use of the techniques provided 
by the Malliavin calculus and, more precisely, we will apply Bouleau-Hirsch's 
criterion (see [I] or Theorem 2.1.2 in [5]). 

First of all, we describe the Gaussian context in which we will use the 
tools of the Malliavin calculus. Namely, we consider the Hilbert space Ht = 
L'^{[0,T];7i) and the Gaussian family of random variables {W{h), h G Ht) de- 
fined at the very beginning of Section [3] Then (W{h),h G Ht) is a centered 
Gaussian process such that E{W{hi)W{h2)) = (/ii, /i2)-Ht) hi,h2 € Ht, and 
we can use the differential Malliavin calculus based on it (see, for instance, [9]). 
The Malliavin derivative is denoted by D and, for any > 1, the domain of 
the iterated derivative in LP(rj;Hf^) is denoted by B^^p, for any p > 2. 
We shall also use the notation 



The first step in order to apply Bouleau-Hirsch's criterion is to study the 
Malliavin differentiability of u{t, x), for all fixed {t, x) £ (0, T] x R^. Recall that, 
for any random variable X in the domain of the derivative operator D, DX 
defines an Ti^— valued random variable. In particular, for some fixed r e [0,T], 
DX{r) is an element of H, which will be denoted by DrX. In the sequel we will 
use the notation ■ and * to denote, respectively, the time and H variables. 




and for d > 3, (I4.13P holds if and only if 




= np>i nfc>i B'^-p. 
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Proposition 5.1 Assume that T satisfies Hypothesis A. Suppose also that the 
coefficients b and a are functions with bounded Lipschitz continuous deriva- 
tives. Then, for any {t,x) G [0,T] x R'', u{t,x) belongs to 'D^'P, for any 
V G [l,oo). 

Moreover, the Malliavin derivative Du{t,x) defines an TIt — valued process 
that satisfies the following linear stochastic differential equation: 

Dru{t, x) —a{u{r, *))T{t — r,x — *) 

+ 11 T{t~ s,x-y)a'{u{s,y))Dru{s,y)W{ds,dy) 

Jr JRt' 

+ [ [ b'{u{s,x - y))Dru{s,x - y)T{t - s,dy)ds, (5.14) 

for all r £ [0, T] . 

The stochastic integral on the right hand-side of Equation (|5.14p must be 
understood by means of the Hilbert-valued integration setting described at the 
very final part of Section [3l 

Concerning the Hilbert-valued pathwise integral, it is defined as follows: 
let ^ be a Hilbert space, {ej)j>i a complete orthonormal system of A and 
{Y{s,y), {s,y) £ [0,T] x U.'^} an .A— valued stochastic process such that 

sup < -foo. 

(s,j/)e[0,T]xR'' 

Then, the A— valued integral 

It= f f Y{s,y)T{s,dy)ds 

Jo JR'' 

is determined by the components ^ J^^ {Y{s, y), ej)^^ r(s, dy)ds, j > ij , which 
are real-valued integrals. Moreover, one can obtain an upper bound for the 
moments of the above integral (see [H], p. 24): 

E{\It\P) < f sup E{\\Y{s,z)rj, [ T{s,dy)ds, p>2. (5.15) 

Eventually, notice that owing to Proposition 13.31 the first term on the right 
hand-side of (|5.14p is fully defined. 

Proof of Proposition \ 5.1\ The statement is almost an immediate consequence 
of Theorem 2 in [T3]. Indeed, the authors of this latter reference prove that, 
under the standing hypothesis and for any fixed {t, x) G [0, T] x W^, the random 
variable u{t,x) belongs to B^'^, for any p € [l,oo). In addition, they show 
that there exists an TYt— valued stochastic process {Q(t, x), {t, x) G [0, T] x R''} 
satisfying 

sup E{\\eit,x)rn^) <^ 

(t,x)e[0,T]xm.'^ 
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and such that, in TIt, 

Du{t,x) ^Q{t,x) + [ [ r{t-s,x-y)a'{u{s,y))Du{s,y)W{ds,dy) 
Jo Jr'' 

+ / b'{u{s,x — y))Du{s,x — y)r{t — s,dy)ds. 
Jo Jr'' 

Moreover, it holds that 

E{\\Qit,x)\\j^^) = *M"(-, 

Hence, in order to conclude the proof, we only need to show that the Hilbert- 
valued random variables Q{t, x) and T{t — •, a; — *)(t(u(-, *)) coincide as elements 
of Ht- 

Let (r„)„>i be the family of smooth functions defined in the proof of Lemma 
13.21 Then, in the proof of Theorem 2 in the process 8 is defined by the 
following limit in TCt'. 

Q{t, x) =Ht - lim Tn{t - ■,x - *)(t(u„(-, *)), 

n — ^oo 

where {u,,(t,x), {t,x) e [0, T] x M''} is the unique mild solution to an equation 
of the form (|1.3p but replacing T by r„. 

As a consequence of the proof of Proposition 3 from |T3] , it is readily checked 
that 

hm ii;(||r„(t ~;x- *)[aK(., *)) - a{u{; = 

and 

\im E{\\[Tr.{t -■,x-*)-T{t-;x- *)]a{u{; *))||?, J = 0. 

Thus, we get that Q{t, x) = a{u{t — •, a: — *))r(t — •, a; — □ 

The main result of the section is the following: 

Theorem 5.2 Assume that T satisfies Hypothesis A. Suppose also that the co- 
efficients b and a are functions with bounded Lipschitz continuous derivatives 
and that |cr(2:)| > c > 0, for all z £ and some positive constant c. Then, for 
all t > and x G M'', the random variable u{t, x) has an absolutely continuous 
law with respect to Lebesgue measure on M. 

Proof. Owing to Bouleau-Hirsch's criterion and Proposition 15.11 it suffices to 
show that \\Du{t,x)\\-HT > almost surely. 

To begin with, from Equation (|5.14p we obtain 

/ \\Dsu{t,x)f^ds>l f \\T{t-s,x-*Mu{s,*mi^ds-I{t,x;S), (5.16) 
Jo ^ Jt-S 



14 



for any S > sufficiently small, where 



t-s 



T(t — r,x ~ z)a'{u{r, z))Dsu{r, z)W{dr, dz) 



r(t — r, dz)h' {u{r, x — z))Dsu{r, x — z)dr 



ds. (5.17) 



H 



The above term /(t, x] 5) may be bounded by 2(/i(i, x; 5) + /2(ij 2;; 5), with 



l2{t,x;S) 



T{t — r,x ~ z)a'{u{r, z))Dt-su{r, z)W{dr, dz) 



T{t — r, dz)h' {u{r, x — z))Dt_gu{r, x — z)dr 



H 



ds, 

H 

(5.18) 
ds. 

(5.19) 



In order to bound from below the term in the left hand-side of (|5.16p . let us first 
obtain a lower bound for the first one on the right hand-side. For this, we will 
make use of the family of smooth functions (r„)„ and ( J^'^)n, considered in the 
proof of Lemma [3?2l that regularise the measures F and F(-, x — *)a{u{t — ■, *)), 
respectively. Then, by the proof of Lemma 13.21 the very definition of the norm 
in Tis and the non-degeneracy assumption on a, we have 



t-5 



\\T{t -s,x- *)a{uis,*))\\i^ds = ||F(.,x - *)a(u(i - •, 



n — 'oc 



lim 

n — *oo 







Jn^'is, y)f{y - z)Jl;''{s, z)dydzds 



> c lim / / / Tn{s,x - y)f{y ~ z)T.a{s,x - z) 
= hm ||F„(.,x-*)||2 =c2||F(.,x-*)||2 =c25(,5), 



(5.20) 



where 



j{S) / / \:Fr{sm\'t^{dOds. 

Jo Jw 



(5.21) 



Now we find out upper bounds for the expectation of the terms Ii{t,x;d) 
and l2(t,x]S). First, in order to deal with the former term, one can use the 
bound (|3.1ip . Thus, taking into account that a' is bounded, we get the following 
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estimate: 

E{h{t,x-5)) 

<C sup i?(||A-.w(t-r,y)||2,J / / \TV{sm\^ti{dOds 

(T,a)e(0,(5)xR'' Jo JR<* 

sup E{\\Dt-.u{t-T,y)\\l^;)g{5). (5.22) 

(T,a)e(0,(5)xR'' 

On the other hand, by (|5.15p the term E(l2{t, x; 5)) corresponding to the Hilbert- 
valued pathwise integral can be bounded by 

Eih{t,x;S))<C sup E{\\Dt-.u{t-r,yWn,)h{S), (5.23) 

(r,a)e(0,5)xR<i 

where h{S) = J^^ r(s, dy)ds. 

At this point, we wiU make use of the foUowing fact (see Lemma 5 in [14j 
and [12], p. 53): 



sup 



E[\\Dt^.u{t-r,y)f^^) <C{g{S)y, (5.24) 



for any g> 1. Hence, bv ^^5^22^ and ([QS]) the terms E(L it. x: 5)) and E(h(t. x: 5)) 
may be bounded, up to constants, respectively, by ((7((5))^ and g{6)h{5), which 
implies that 

E{I{t, X- 5)) < CgiS)ig{S) + HS)). (5.25) 
For any fixed small 6 > 0, let n be a sufficiently large positive integer such 
that i < ^g{S). Then, owing to (|5.20p and (|5.25p and applying Chebyshev's 
inequality, we obtain 

P (^^ p,zi(i,x)l|2,rf,s < ^) < ^ (^I{t,x;6) > jg{6) - 

< (^gi^)--) E{I{t,x-6)) 



<(^y5(5)--j gmg{5) + h{5)). (5.26) 
Therefore 

\\m p( [ \\DM-t. x)\\nds < < ^(^W + K^)): 
and the latter term converges to zero as 5 tends to zero. Hence, 
p(^j^ \\Dsu{t,x)\\l^ds = {)^ =0, 

which concludes the proof. □ 

Remark 5.3 In the particular case of the three-dimensional stochastic wave 
equation, the above Theorem \5.2\ generalises Theorem 3 in the reference US]j . 
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6 Smoothness of the density 

This section is devoted to prove that, for any fixed {t,x) e (0,T] x M'', the 
law of the random variable u{t,x) has an infinitely differentiable density with 
respect to Lebesgue measure on R. This will be achieved by showing that u{t, x) 
belongs to the space D°° and that the inverse of the Malliavin matrix of u{t, x) 
has moments of all order (see, for instance. Theorem 2.1.4 in 9 ). 

Recall that for any differentiable random variable X and any iV > 1, the 
iterated Malliavin derivative X defines an element of the Hilbert space 
As for the case TV = 1, for any r = (ri, . . . , tat) € [0,r]^, the 
element DX[r) of W^^ wiU be denoted by DrX. We will also use the notation 

for e [0, T], (y9i G i = 1, . . . , TV. In particular, we have that 

dri...drN l^((n,e,J,...,(riv,e,„))^P, 

where {ej)j>o is a complete orthonormal system of Ti.. Let 

A^{g,X) ■.^D^g{X)^g'{X)D^X, 

where a = ((ri, yji), . . . , (rjv, (/jat)), £ [0,T] and ipi G H. Notice that 
{g, X) = if = 1 and it only depends on the Malliavin derivatives up 
to the order TV - 1 if > 1. 

We now state the main result concerning the Malliavin regularity of the 
solution u{t, x). 

Proposition 6.1 Assume that T satisfies Hypothesis A. Suppose also that the 
coefficients a and b are C°° functions with bounded derivatives of any order 
greater than or equal to one. Then, for every {t, x) G [0, T] X M'', the random 
variable u{t,x) belongs to the space D°°. 

The iterated Malliavin derivative D^u{t,x) satisfies the following equation 
in LP{n; H^^), for any p > 1 and N > 1: 

D^u{t,x) = Z'^{t,x) 

+ [ [ T{t- s,x - z)[S^{a,u{s,z)) + D^u{s,z)a'{u{s,z))]Wids,dz) 
Jo Js.'^ 

+ [ ds [ r{s,dz)[S^{b,u{t~s,x-z)) 
Jo Js.'' 

+ D^u{t -s,x- z)b'{u{t - s,x~ z))], (6.27) 
where Z^{t,x) is the element of LP{^l;'H®^) defined by 

N 

(Z^{t,x),ej^ ® • ■ • (g) ej^)-H<sN = ^(r(t -r,,x- *)D^-^a{u{r,, ej^)H, 

4=1 

(6.28) 



vl|2 
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for any r = (ri, . . . ,rAr) G [0,r]^ and ji, . . .Jn G {1, ■ • ■ ,^}- 

Moreover, it holds that 

sup E{\\D^u{s,y)\\l 

«n) < +00, 

for all p > 1 . 

Proof. It is almost an immediate consequence of Theorem 1 in [14| and Propo- 
sition [O] from the preceding Section [5l 

Namely, we just need to check that, using the same notation as in the proof 
of Theorem 1 in [T3], the sequence of L^(17; 7Y®^)— valued random variables 
(Z^'"(i,a;))„>i converges to Z^{t, x) as n tends to infinity. We should mention 
that in that reference, Z^''"-{t, x) is constructed using a regularisation procedure, 
that is smoothing the measure F, and by means of a similar expression to (|6.28p . 
In [14], the objective was to define the initial condition of the linear stochastic 
equation satisfied by the iterated Malliavin derivative D^u{t, x) as the limit of 
Z^''^{t,x). We claim that this limit equals to Z^{t,x), defined in the present 
proposition's statement (see (|6.28I) V 

Indeed, the convergence of Z^'"(t,x) to Z^(t,x) in L^{n;n^'^) can be 
easily studied using the same arguments as in the proof of Lemma 3 in |14| . □ 

We are now in position to state and prove the main result of the paper. 

Theorem 6.2 Assume that T satisfies Hypothesis A, the coefficients a and b 
are C°° functions with bounded derivatives of any order greater than or equal to 
one and \cr{z)\ > c > 0, for all z G R. Moreover, suppose that there exist 7 > 
such that for all t G (0, 1], 

f I \TY{sm\^y.{di)ds>C^T\ (6.29) 

for some positive constant Ci. Then, for all {t,x) G [0,r] x M.'^, the law of 
u{t, x) has a C°° density with respect to Lebesgue measure on R. 

Proof. In view of Proposition 16.11 we need to show that the inverse of the 
Malliavin matrix of u{t, x) has moments of all order, that is 




for all q > 2. It turns out (see, for instance, Lemma 2.3.1 in [9J) that it suffices 
to check that for any q > 2, there exists an £o(<z) > such that for all e < eo 




(6.30) 
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Proceeding as in the proof of Theorem 15.21 for any S > sufficiently small we 
obtain the following estimate: 



P 



xWnds <e] <P{ I{t, x; S) > -g^S) - e 



g{5)-e] E{\I{t,x-,5)\P), (6.31) 



< 



for any p > 0, where we recall that /(t, x; 5) is defined by (|5.17p and g{5) is 
given by ((OT|) . 

We decompose now the term /(i, x] 5) as in the proof of Theorem 15.21 so 
that we need to find upper bounds for E{\Ii{t,x;5)\Y , i = 1,2 (see (|5.18l) and 
(|5.19p ). On one hand, owing to Holder's inequality and p. lip we get 



E{\h{t,x;S)n 



= E 



< SP-'E 




T{t — r,x ~ z)a'{u{r, z))Dt-su{r, z)W{dr, dz) 



ds 



H 



T{t — r,x — z)a'{u{r, z))Dt-su{r, z)W{dr, dz) 



2p 

ds 

H 

(6.32) 



<SP-\gi6)r sup E(\\D,^.uit-r,y)\\X 

(t,j/)6[0,T]xR'' ^ 

The above estimate ()6.32p let us conclude that 

Ei\h{t,x;S)n<CSP-\9{S)r. 

On the other hand, using similar arguments but for the Hilbert-valued pathwise 
integral (see (|5.15l) '). one proves that E{\l2{t,x;6)\P) may be bounded, up to 
some positive constant, by SP~^{g{6)y. Thus, we have proved that 



P 



(^J^ \\Dsu{t, x)r^ds <e^<C (^^g{5) - ' SP-\g{S)r. (6.33) 



At this point, we choose d — (5(e) in such a way that g{S) 
implies that 4-e > CS'^, that is (5 < Ce~ . Hence, 



\\Dsu{t,x)\\l^ds <e] < Ce~ 



By (jOO)) . this 



and it suffices to take p sufficiently large such that > 



□ 



Remark 6.3 As it is pointed out in USf . Appendix A (see also '61), when T 
is the fundamental solution of the wave equation in M.'^, with d = 1,2,3, then 
condition \6.29j) is satisfied with 7 = 3. 

On the other hand, if V is the fundamental solution of the heat equation on 
, d>l, then condition i6. 29]) is satisfied with any 7 > 1 (see Lemma 3.1 in 
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Remark 6.4 Theorem \ 6.2\ vrovides a generalisation of Theorem 3 in JI^ for 
the case of the three-dimensional wave equation (see also |j5j/). Moreover, it 
also generalises the results in |^ for the stochastic wave equation with space 
dimension d = 1,2 and the stochastic heat equation in any space dimension. 
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